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Abstract. Proper vertex colorings of a graph are related to its boundary map d\, also called its signed 
vertex-edge incidence matrix. The vertex Laplacian of a graph, L = did[ = D — A, a. natural extension of the 
boundary map, leads us to introduce nowhere-harmonic colorings and analogues of the chromatic polynomial 
and Stanley's theorem relating negative evaluations of the chromatic polynomial to acyclic orientations. Further, 
we discuss several examples demonstrating that nowhere-harmonic colorings are more complicated from an 
enumerative perspective than proper colorings. 



1. Introduction 

A classical line of investigation in graph theory is the study of proper graph colorings. Two well-known 
enumerative theorems regarding such colorings are the Birkhoff-Whitney theorem establishing the existence 
of the chromatic polynomial and Stanley's reciprocity theorem, given as Theorems 2. 1 and 2.2 below. Proper 
vertex colorings of a graph are related to its boundary map di , also called its signed vertex-edge incidence 
matrix; this relationship was exploited by the first author and T. Zaslavsky in [2] to produce geometric proofs 
of Theorems 2. 1 and 2.2. In this paper, the vertex Laplacian of a graph, L = d[d[ = D —A, a natural extension 
of the boundary map, leads us to introduce nowhere-harmonic colorings. We prove Theorems 3.4 and 3.5, 
intriguing analogues of Theorems 2.1 and 2.2. Further, we discuss several examples demonstrating that 
nowhere-harmonic colorings are more complicated from an enumerative perspective than proper colorings. 

Our paper is organized as follows. In Section 2 we establish notation and review proper graph colorings, 
including Theorems 2.1 and 2.2. In Section 3 we introduce nowhere-harmonic colorings and state our main 
results. In Section 4 we discuss the theory of inside-out polytopes. In Section 5 we provide proofs of 
our theorems. In Section 6 we provide the results of some computer calculations and discuss in detail the 
nowhere-harmonic enumerators for star graphs. 

2. Proper Graph Colorings 

Let V = ,X2, . . . ,x„} be an ordered set, with the order on V coinciding with the order on the subscripts. 
A graph G on V is a pair iy,E) where £" is a collection of 2-element subsets of V . As usual, we call the 
elements of V vertices of G and the elements of E edges of G. An m-coloring of G is a labeling of the 
vertices of G with values from {1, . . . ,m} and a proper m-coloring of G is a labeling of the vertices of G 
with labels from the set {1, ... ,m} such that adjacent vertices have distinct labels. We denote by Xg{^i) the 
number of proper m-colorings of G. The following theorem is the first prerequisite for our investigations. 

Theorem 2.1 (Birkhoff [5], Whitney [16]). For every graph G with n vertices, Xg{^) given by a polyno- 
mial in m of degree n. 
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We call the polynomial Xcifn) the chromatic polynomial of G. A priori, the value of this polynomial has 
a combinatorial interpretation only for positive values of m; a beautiful theorem providing a combinatorial 
interpretation for negative integral values of m is the second prerequisite for our investigations. 

We refer to the (edge) orientation of G given by directing edge {x,-,xy} as {xt < xj} whenever / < j as 
the standard orientation of G. We say an arbitrary orientation £ of the edges of G is acyclic if it does not 
contain a directed cycle in G. Note that the standard orientation is acyclic. We may represent any orientation 
e as a sign function on the edges of G, i.e., e{e) = 1 if the orientation of e under e agrees with the standard 
orientation and £{e) = —I otherwise. We say a vertex coloring c of G is compatible with e if for every 
directed edge e = {x, < xj} G £ in the standard orientation, either e(e) and c{xj) — c(x/) have the same sign 
or c{xj) — c{xi) = 0. 

Theorem 2.2 (Stanley [14]). For every graph G with n vertices and every integer m > 0, 

(-ir;tG(-m) = £a(c), 

c 

where the sum ranges over all m-colorings c of G and a{c) is the number of acyclic orientations of G 
compatible with c. In particular, {~1)"Zg{—^) equals the number of acyclic orientations ofG. 

One proof of Theorem 2.2 is via the theory of inside-out polytopes developed by the first author and T. 
Zaslavsky in [2]. The theory of inside-out polytopes is a framework for enumerating the points of Z" con- 
tained in the relative interior of a convex polytope ^ c M" but not in a given affine hyperplane arrangement 
J^, where both the polytope and hyperplane arrangement are rational. One may view inside-out theory as 
a template for producing geometric reciprocity theorems for which a combinatorial interpretation may be 
sought. In the case of proper graph colorings, geometry enters the picture in the following way. 

Given a graph G, we denote by Co the real vector space with basis given by the vertices of G; we denote 
by C\ the real vector space with basis given by the edges of G. We may then think of a vertex m-coloring of 
G as an integer point in Cq. It is immediately clear that a proper coloring of a graph G is an integer point in 
Co \ (G), where Jf{G) is the hyperplane arrangement induced by G, namely 

(G) := [xi =Xj : {xi,Xj} G £} , 

where a hyperplane arrangement is a finite collection of affine hyperplanes in some Euclidean space. This 
setup is the starting point of [2] for studying the chromatic polynomial and acyclic orientations of G from a 
geometric point of view; we will formulate this in a more general way as follows. 

For every graph G= iy,E), there is a natural linear boundary map 5i : Ci ^ Co, where di {e) = xj — x; if 
e = {xi < Xy } G £. If we represent d\ as a matrix, then its transpose d[ : Co ^ Ci is the coboundary map. We 
will liberally use the same notation for a linear transformation and its matrix representation. If / : R" — > is 
a linear transformation, then a hyperplane arrangement J>^ in M*^ induces a hyperplane arrangement (/) 
in M", the pullback of Jif, as follows: any hyperplane /i G is given by a linear functional (p on M!^, and 
thus <p o/ is a linear functional which, in turn, defines a hyperplane in W; the collection of these hyperplanes 
is "(/). (Fullbacks were used in [2] to study antimagic labellings.) We denote by = {xe = 0:e EE} 
the arrangement of the coordinate hyperplanes in Ci (or in Co; which space is being referenced will be clear 
from context). 

Lemma 2.3. The proper colorings ofG are precisely the positive integer points in Cq \ 

Proof. We need to show that JfQ{d[) = Jf{G). The hyperplane ft : Xg = in induces the hyperplane 
d{{<p) = in J^{d{), where <p is the linear form associated with h. But 5j o <p is the zero matrix except for 
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the row, which is identical to the row of d[. Thus the hyperplane given by ^{(<jp) = is xj —xi = 0, 
where e = {xi ,Xj}. □ 

Thus, proper colorings of a graph G arise naturally as positive integer points in Co that avoid the hy- 
perplane arrangement J^^{d{). While the combinatorial and geometric importance of the boundary map 
for G is readily apparent, the relation between the coordinate arrangement J^q in C\ and G is less appar- 
ent. The importance of the coordinate arrangement can be seen by noting that there is a natural one-to-one 
correspondence between regions of and orientations of the edges of G, via the following construction: 

• Mark each region of J^q with a sign vector s G {±1}'", where each sign indicates whether this 
particular coordinate is larger or smaller than 0; 

• Apply s to the standard orientation of the edges of G, negating an edge precisely if the corresponding 
coordinate of s is negative. 

This construction can be reversed and thus gives the desired correspondence. (This idea is not new and goes 
back to Greene.) As demonstrated in [2], it is this correspondence that is at the heart of Theorem 2.2. 

In addition to proper colorings, nowhere-zero flows and tensions on graphs arise from a geometric per- 
spective by studying the boundary map and coordinate hyperplane arrangements. Enumerative properties of 
these structures have been the subject of previous papers, such as [3, 6, 7, 9, 10]. 



3. Nowhere-Harmonic Colorings and Main Results 

When in an enumerative situation where inside-out theory applies, one automatically obtains a geometric 
reciprocity statement that parallels Theorem 2.2. However, inside-out theory does not immediately yield a 
combinatorial interpretation for this reciprocity statement (e.g. [2, Problem 7.2]). One of the contributions 
of this work is to demonstrate that in the case of nowhere-harmonic colorings of graphs, we can find a 
such a combinatorial interpretation via Theorem 3.5. We will begin by providing an arithmetic definition of 
nowhere-harmonic colorings. 

Definition 3.1. An m-coloring c of G is nowhere harmonic if for every vertex v of G, the value of c(v) is 
not equal to the average value of c(w) where the average is taken over all w that are neighbors of v. We 
denote by hcim) the number of nowhere-harmonic m-colorings of G. We denote the generating function for 
nowhere-harmonic colorings as 

Hg{z) := £ hG{m)f. 

m>\ 

Two basic properties of nowhere harmonic functions and their enumerators are the following. 

Proposition 3.2. Every graph G admits a nowhere-harmonic 2-coloring. 

Proposition 3.3. IfG has at least one edge, then Hg{z) is divisible by 2z}. 

Given Theorem 2.1, it is natural to ask if the function hcim) has nice structure. One answer is the 
following; recall first that a quasipolynomial is a function /(m) defined on Z such that there exist A: G N and 
polynomials pi{m),p2{m), . . . ,Pk{m) such that /(m) = Pj{m) whenever m = j mod k. The minimal such k 
is the period of /(m) and the polynomials pi{m),p2{m), . . . ,Pk{m) are its constituents. 

Theorem 3.4. For every connected graph G with n vertices, there exists a quasipolynomial f(m) of degree 
n such that hcim) = f{m). Further, hcim) has ( constant) leading coefficient 1 and |^g(~ 1 ) I = 2" — 2. 
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Due to the fact that hcim) is a quasipolynomial, Hg{z) is a rational function of the form 

gjz) 

(l_^p)n+l' 

where p is the period of %Q{m); for justification see [15, Chapter 4] or [1, Lemma 3.24]. As our examples 
in Section 6 indicate, there are subtle and interesting cancellations that occur in these rational forms. 

To state our analogue of Theorem 2.2, we need a few preliminary definitions. Given a graph G, a coloring 
c of G, and a vertex v of G of degree dy, we say c is subharmonic at v if 

c(v) < ^ 

{v,w}^E{G) 

and we say c is superharmonic at v if 

c(v) > J- £ c(w). 

{v,w}e£(G) 

We call a it 1 -labeling of the vertices of G a vertex orientation of G and say that the vertex orientation is 
proper if it is not the all-ones labeling or the all-negative-ones labeling. Note that since G has n vertices and 
two choices of sign for each vertex, there are 2" — 2 proper vertex orientations of G. Let £ be a vertex orien- 
tation of G. Given a coloring c of the vertices of G, we say that c is compatible with e if c is subharmonic at 
V when e(v) = 1 and c is superharmonic at v when e(v) = — 1. 

Theorem 3.5. For every connected graph G with n vertices and every integer m> 0, 

(-iyho{-m)=l^P{c), 

c 

where the sum ranges over all m-colorings cofG and j3(c) is the number of proper vertex orientations ofG 
compatible with c. 

As with Theorem 2.2, Theorem 3.5 is a consequence of the theory of inside-out polytopes where the role 
of the boundary map of a graph is played by the vertex Laplacian of a graph. 

Definition 3.6. Let G = (V,^) be a simple, loopless graph with n vertices. Let the elements of V be ordered 

asV = {xi,X2, ...,Xn}. The vertex Laplacian (also called the Kirchhojf matrix) of G is 

L:=D-A, 

where Z) is the n x n matrix indexed by the elements of V such that Z), y is the degree of if i = j and 
otherwise and A is the 0/1-valued adjacency matrix of G indexed in the same way as D. 

The following theorem connects the vertex Laplacian to the boundary map. 

Theorem 3.7. (See [4, Chapter 6].) If di is the boundary map of G, then 

(1) L=d,dl 

hence L is a positive semi-definite real symmetric operator and the dimension of the nullspace ofL is equal 
to the number of connected components ofG. 
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While the following is well known to algebraic graph theorists, it is helpful for others to compare the 
vertex Laplacian for a graph G with the usual Laplacian for domains in the complex plane, A := ^ + 
This is an appropriate analogy, as we can write 



div(grad) 



d d 
3^ 



d 
dy 



which corresponds to the factorization (1). Thus, the transpose of the boundary map for a graph is a discrete 
analogue of the gradient operation for functions on domains while the boundary map for a graph is a discrete 
analogue of the divergence of a vector field. 

A smooth function / on a domain in the plane that satisfies A/ = is called a harmonic function. Given 
a graph G with vertex Laplacian L, it is not hard to show that if c is a harmonic function on G, i.e., if Lc = 0, 
then c is constant on the connected components of G and that if G is connected, then the nuUspace of L is 
spanned by the all-ones vector. More generally, one has the following combinatorial interpretation of the 
equation Lc = b: if we consider c as a labeling of the vertices of G, then the entries in b are the differences 
between c{xi) and the average of the values of c{xj) for all j such that Xj is adjacent to Xi. Thus, we have the 
following lemma. 

Lemma 3.8. An m-coloring of a graph G is nowhere-harmonic if and only if it corresponds to a positive 
integer point in Co \ (L) = Cq \ J^q {d\ d[). 

Given the intimate connection between the boundary map and proper graph colorings. Theorem 3.7 indi- 
cates that the vertex Laplacian is a natural next object to consider from an enumerative point of view, hence 
our study of nowhere-harmonic colorings. It should be noted that while harmonic functions on finite graphs 
are simply described, harmonic functions on infinite graphs have a rich structure and are the subject of active 
research, see for example [13] and [12]. 

The importance of the coordinate arrangement to nowhere-harmonic colorings can be seen by noting that 
there is a natural one-to-one correspondence between regions of J% in Co and =bl-labelings of the vertices 
of G, via the following construction: 

• Mark each region with a sign vector 5 G {±1}'", where each sign indicates whether this particular 
coordinate is larger or smaller than 0; 

• The resulting vectors s obviously correspond bijectively to itl-labelings of the vertices of G. 

By specifying a 1 or —1 at each vertex, we are specifying whether that vertex should be a host for sub- or 
superharmonic behavior for colorings, much as the orientation of an edge restricts the labeling a coloring 
may have in Theorem 2.2. 

Remark 3.9. A nowhere-harmonic coloring is not necessarily proper. It might be interesting to investigate 
enumerative properties of colorings that are both proper and nowhere harmonic. 

4. INSIDE-OUT POLYTOPES 

Our proofs of our theorems will involve inside-out theory. Given a convex polytope ^ cM." and a 
hyperplane arrangement Jf , where both the polytope and hyperplane arrangement are rational, we call 

{I^,Jif) an inside-out polytope. A region of Jif is a connected component of W \ J^. A closed region is 
the closure of a region. A region of M') is the nonempty intersection of a region of M' with A vertex 
of {^,J^) is a vertex of any such region. A closed region of {^,J€^) is the closure of an open region of 
M') and therefore meets the relative interior The denominator of M') is the minimum pel>>i 
such that the vertices of the p-ih dilate of {^,Jif) all have integer coordinates. 
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Given a point v G W, the multiplicity of v with respect to (,^, ), denoted m^^jf (v), is the number of 
closed regions of {^,J^) containing v if v G <^ and otherwise. The closed Ehrhart quasipolynomial for 
(■^, J^) is the function defined for ? G Z>i given by 

ve(f-iZ")n.^ 

The open Ehrhart quasipolynomial for ) is the function defined for ? G Z>i given by 

£^.,^(0:=#(riz"n[^°\jr]). 

The main result of inside-out theory is the following. 

Theorem 4.1, [2, Theorem 4.1] If is a closed, full-dimensional, rational inside-out polytope in 

M", then Eff>^jfe{t) and £"^0^(0 are quasipolynomials in t that satisfy the reciprocity law E'^^ jf{t) = 
{—l)"Egg^jjf{—t), with period equal to a divisor of the denominator of{^, J^f), with leading term vol(=^)?", 
and with constant term E^^jfiO) equal to the number of regions of{^,Jf). 

The main idea behind the proof of Theorem 4. 1 is that the arrangement dissects into disjoint, open, 
rational polytopes. Ehrhart theory, see [1] or [15, Chapter 4], then implies that the function counting lattice 
points in dilates of a given region /? is a quasipolynomial, called the Ehrhart quasipolynomial ofR, whose 
degree is equal to the dimension of the region and whose period divides the denomenator of R. By summing 
up the Ehrhart quasipolynomials of each open region, we recover the counting function E^o _^{t)- One 
may also compute the function E°^o j^{t) via Mobius inversion on the intersection lattice for the hyperplane 
arrangement. The reciprocity law in Theorem 4. 1 is then a consequence of Ehrhart-Macdonald reciprocity 
for rational polytopes. 

5. Proofs 

Proof of Proposition 3.2. We may assume that \V\ > 2. With two colors, the only way that a coloring could 
avoid being nowhere harmonic would be for some vertex to have the same label as all of its neighbors. Also, 
note that we only need show that we can produce such a coloring for connected graphs, as disconnected 
graphs can be dealt with component by component. Let T be a rooted spanning tree of G with root t. Label t 
with 1. Label all the neighbors of t with 2. Label all of their neighbors with a 1. Proceed inductively until all 
vertices have been exhausted. When this labeling is complete, every vertex v of G has at least one neighbor 
with a different label from v, hence our coloring is nowhere harmonic. □ 

Proof of Proposition 3.3. Observe that for any graph G with n vertices and Laplacian L, we have hoim — 
1) = E°^ jr''(L)^"^-'' ^ denote the inside-out polytope ((0, 1)", J^''(L)). Let T be the invertible affine 
transformation given by x —x + E, where E denotes the vector with all entries equal to 1. Since —x 
is in the cube [—1,0]" and E is in the kernel of L, T takes £P onto itself. As T is invertible, the bijective 
correspondence between the regions follows. Thus, every nowhere-harmonic m-coloring c of G has a natural 
paired coloring T{c), and hence Hciz) is divisible by 2. As G has at least one edge, G does not admit 
nowhere-harmonic 1-colorings; thus, the first term in Hg{z) is the term and our claim is verified. □ 

Proof of Theorem 3.4. Given that for any graph G with n vertices and Laplacian L, we have hcim — I) = 
E° J (m). Theorem 4.1 immediately imphes the quasipolynomiaUty and degree results for hoim). To 

(0,l)",Jt^ (L) 

show that |/ig(-1)| =2"-2, we need to prove that the number of regions of ([0, 1]",^''(L)) is2"-2. This 
is a straightforward exercise in hnear algebra, as follows. 
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As L is a positive semi-definite real symmetric operator, the spectral theorem implies that L has only 
noimegative real eigenvalues and an orthogonal basis of eigenvectors. Since G is connected, we can conclude 
that the all-ones vector spans the kernel of L, and that L is invertible on Ker(L)-'-. Thus, each hyperplane K 
in the arrangement .J^q{L) is spanned by Ker(L) and fi Ker(L)-'-. The intersection of K with Ker(L)-'- is 
the inverse image of one of the coordinate hyperplanes in Co restricted to lm(L). As L restricted to Ker(L)-'- 
is invertible with positive eigenvalues and orthogonal eigenvectors, we can conclude that the regions in 
([0, \Y ,^q{L)) are in one-to-one correspondence with the regions in Im(L) n J^. There are exactly 2" — 2 
such regions, since Im(L) intersects all the open regions defined by the coordinate arrangement except the 
positive and negative orthants. □ 

Proof of Theorem 3.5. Our goal is to interpret the geometric reciprocity portion of Theorem 4.1 combinato- 
rially in this setting. From the proof of Theorem 3.4, we know that 

and that the regions of ( [0 , 1 ] " , J^" (L) ) correspond bij ectively with proper vertex orientations of G. To finish 

the proof, note that E, ,„ y/--,7\i>^)y by definition, sums over the colorings described in the theorem and the 

[U.IJ .-Xg (L) 

multiphcity of each coloring c in the sum is exactly j3 (c). □ 



6. Examples 

Using the I DP package [8], we computed sample generating functions for nowhere-harmonic quasipoly- 
nomials among four classes of graphs: path graphs, cycle graphs, complete graphs, and star graphs. In 
addition to our computations, we provide in this section a detailed analysis of nowhere-harmonic colorings 
of star graphs to show how polyhedral geometry can illuminate the structure of these generating functions. 

To put the following algebraic expressions into context, note that for a rational polytope Q with vertex 

coordinates, indexed by the set /, of the form < — :/£/> one would hope that the period of the Ehrhart 

U« J 

quasipolynomial for Q is lcm{^,- : i El} (which is the maximal possible period, see, e.g., [1, Section 3.7]). 

However, for general rational polytopes, the period of an Ehrhart quasipolynomial can be smaller than the 
1cm of the denominators of the vertex coordinates, so one must exhibit some caution when drawing geomet- 
ric conclusions from the algebraic expressions. Nonetheless, when studying (inside-out) polytopes arising 
from combinatorial constructions, one hopes that the denominators of the relevant vertices are reflected in a 
reasonable way in the denominators of rational generating functions. 

In the examples that follow, though Hg{z) is a rational function of the form -^j^^y+r, where p is the period 
of fiQ{m), this is not necessarily the reduced form of the quotient. We have chosen to represent these rational 
functions in reduced form to emphasize the divisors of p rather than p itself. This will prove particularly 
useful when we analyze nowhere harmonic functions of star graphs. The numerators of these expressions 
remain quite mysterious. 



6.1. Path Graphs. The path graph on k vertices has the edges {j,j +1} for 1 < j < k. We denote its 
generating function for nowhere-harmonic m-colorings by Pk{z). 
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2z^(l+2z + 3z^) 

(l-z2)(l-z)3 



2z^ (2 + lOz + 3 + 40z^ + 36z^ + 18z^ + Tz*^) 

(1-Z)2(l-Z2)2(l_^3) 



P5(Z) = 



r 3 + 30z + 173z2 + 480z^ + 979z^ + 1456z^ + 1740z^ ^ 
^ \ +1586z' + 1 185z^ + 648z^ + 273^^° + 72z" + I5z^^ ^ 

(1-Z)(l-Z2)2(l-^3)2(1_^) '- 



f 5 + 82z + 748^2 + 3326z^ + 10964z^ + 27767z^ + 59103z^ + 106774z'' 



2z'< 



+ 171025z^ + 24277U^ + 3124012^° + 363363z" + 386862^1^ 
+ 374277zi^ + 331925zi^ + 266383zi^ + 194327^^^ + 126327zi^ 



+ 73219z'*^ + 36387z''^ + 15522^2" 



5208z2i + 1372^22 + 231^23 + 31^24 



(l-z2)(l-z3)2(i_^4)2(i_25)(i_^6) 



6.2. Cycles. The cycle on k vertices has the edges {j,j + 1}, for I < j < k, and {1,/:}. We denote its 
generating function for nowhere-harmonic m-colorings by 



C3iz) 



6z^(l+z2) 

(l-z)3(l-z2) 



C4(Z) 



2z2(3 + 1 Iz + 41z2 + 42z^ + 29z^ + 1 Iz^ + 7z^) 

(l-z)2(l-z2)2(l-z3) 



Csiz) 



^ { l+7z + 38z^ + 75z^ + 148zV201z5+265z^ + 256z'^ 

^ \ +271z^ + 195z^ + 152z^° + 75z' ^ + 34z^^ + 7z^^ + 3z^'* ^ 
(l-z)2(l-z2)(l-z3)(l-z4)(l-z6) ' 



C6{Z) 



2zH 



' 10+ 108z + 965zV4022z^ + 12814zV31491z^ +65668z^ + 115971z^ 
+ 182775z^ + 254824z^ + 323136z^° + 369837z^^ + 388256z^^ + 369420z^^ 
+ 322890z^^ + 254413z^^ + 182376z^^ + 1 15872z^'' + 65683z^^ + 31509z^^ 
+ 13018z^° + 4109z^^ + 1034z^2 ^ i68z^3 ^ 3j^4 



(1 -Z2)(l -z4)2(l -^3)2(1 _^5)(i _^6) 



6.3. Complete Graphs. The complete graph on k vertices has the edges for 1 < j < j <k. We denote 
its generating function for nowhere-harmonic m-colorings by Kk{z). 
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(l-z)3(l-z2) 

2z^ (7 + 12z + ITz^ + 29z^ + 7z^) 

(l-z)3(l-z2)(l-^3) 

lOz^ (3 + 7z + 28z2 + 19z3 + 50z^ + 5z^ + 32z^ - + 3z^) 
(l-z)4(l-z3)(l-z4) 



2z2. 

Keiz) = 



31 + 179z + 7202^ + 1469z^ + 2400?^ + 3510z^ + 351 k*^ 

+ 3841z' + 2520z*^ + 2220?"^ + 569z"' + 630z^^ - 31z^^ + 31z^^ 



(l-Zr(l-z3)(l-z4)(l_^5) 

6.4. Star Graphs. The star graph on k vertices is the complete bipartite graph K\±^i, i.e., the edges are 
{{l,j} : 2 < j <k}. We denote its generating function for nowhere-harmonic m-colorings by Sk{z). 



S4iz) 



(1-Z)3(1-Z2) 

2z2(7z^ + 12z^ + 26z3 + Hz^ + 9z + 1) 

(l-z)3(l-z2)(l-z3) 

2z2(15z'^ + 50z^ + 177z^ + 265z^ + 323z^ + 309z^ + 187z^ + 95z^ + 18z + 1) 

(1-Z)3(1-Z2)(1-Z3)(1-Z4) 



2z- 



j 3lz^^ + 180z" +9342^2 ^2 174z" +4039z^° + 5980z'' + 7256z^+^ 
\ 1411 z'' + 6386z^ + 4615z^ + 2629z^ + 1 1 89z^ + 324z2 + 45z + 1 



(l-z)3(l-z2)(l-z3)(l-z4)(l-z5) 

Note that a star graph on k vertices is a tree with k—l leaves and a single vertex of degree k—l. On 
the leaves of Ki^^-i^ to satisfy the nonharmonic condition requires only that the labels on the leaves do not 
equal the label on the center of the star, vertex 1. Therefore, in the case of i, every nowhere-harmonic 
coloring is actually a proper coloring (but not vice versa). K\^ic-i is a simple graph to look at from a nowhere- 
harmonic point of view, in the sense that after removing the first row of the matrix L^, one is left with the 
matrix corresponding to the boundary map for Ki^^^i, which is totally unimodular. Thus, by removing 
the hyperplane corresponding to the first row of L^, we are left with an inside-out polytope whose open 
Ehrhart polynomial is actually a polynomial, as all the vertices of the regions are integral. We will see that 
the addition of this single hyperplane creates regions with large denominators leading to these comphcated, 
though not completely chaotic, generating functions for nowhere-harmonic m-colorings. 

There are a variety of interesting patterns to observe regarding the functions Sk{z). As anticipated by 
Proposition 3.3, their numerator polynomials are of the form 2z}f\{z) for some polynomial fi{z). Also, 
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the leading terms of the fk{z)'^ are of the form 2^~^ — 1, which essentially counts regions in the associated 
inside-out polytope. A mystery to us is that the coefficients of the second highest terms of the fk{z)'^ 
appear to be of the form 3*~^ — 2 ■ 2*~^ + 1 = 25(A:, 3), where S{k, 3) is the StirUng number of the second 
kind [11, sequence R\q\. All the denominators of the generating functions in our examples are of the form 



(1 



,/t-l^ 



'1 — z )(1 — z) , though we have not been able to prove that this will always be the case. 



Our goal in the remainder of this section is to investigate the inside-out poly topes k •={{^:^)^: i^k))^ 
where is the vertex Laplacian for Ki^^-i^ with two purposes. First, through a careful analysis of the re- 
gions of we will be able to produce an alternative expression for the generating function Sk{z). Second, 
our analysis will also yield some insight regarding the form of the denominator of Sk{z). 

Fix A: > 3. Observe that the vertex Laplacian for K\^k-i is a block matrix of the form 



1 



1 (k- 
2 



Lk= 3 



1 





1 



1\ 








k\ -I ••• 1 / 

where the rows and columns of are Unearly indexed by {1,2, ... as shown. Each region Rg of the 
inside-out polytope is determined by a proper vertex orientation e of ^Ti^^-i, with the constraint matrix 
for Rg given by 



-In 

where Lj^ is the matrix obtained by multiplying each row of the Laplacian matrix by a — 1 or 1 in accordance 
with e. For example, with L3 and the vertex orientation {li-)-— l,2i-^l,3i-)-— 1}, the corresponding region 
is 





1 

< 1 
1 





Voy 

We will identify the vertices 1,2,3 with variables xj ,X2,X3, . . ., respectively, when discussing equa- 
tions corresponding to rows of the constraint matrix for regions of Figure 1 contains a picture of ^3, 
where the axes are labeled 1,2,3 in accordance with the vertices. Two features of ^3 are immediately ap- 
parent: first, there is a transformation taking the half of =^3 "above" the hyperplane 2xi -I-X2 -I-X3 = to the 
half of ,^3 "below" the hyperplane; this is the transformation T from Proposition 3.3. Second, there are 
only two types of regions, one being pyramids over squares. We can formulate precise analogues of these 
features for general 

We will from now on consider only regions of corresponding to vertex orientations of A'l ^t-i where 
1 I-)- — 1 . We will call these regions negative regions of Our next goal is to produce an alternative 





-2 
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1 \ 






-1 
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-1 
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\ X3 




-1 















-1 
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-1 / 
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Figure 1. The inside-out polytope .^3. 

formula for Sk{z) by utilizing symmetries of Ki^k-\- The automorphism group of Ki^^^i consists of <Sk-i, 
the permutations of the A: — 1 leaves of the star; &k-i induces a permutation transformation of M/^ under 
which has its negative regions permuted. Two negative regions are in the same orbit of this group action 
if they correspond to vertex orientations with the same number of — I's. It follows that there are exactly 
k—l orbits of this ©yt-i -action and, for 1 < j < ^ — 1, the orbit is represented by the region defined 
by negating the first j rows of L^, i.e., the region for the vertex orientation {1 — 1,2i-^ i-> 
— 1 H-> 1, . . . h-)^ 1}. We may now compute ^^(z) by summing the open Ehrhart series for the ^I's 
weighted by the number of elements in each respective orbit, obtaining the following proposition. 

Proposition 6.1. Let be defined as just described, and let 

Siiz) :=f:#(m-iZ«n^/)z'". 

m=l 

The following equality holds: 

Unfortunately, in this case computing Sk{z) is not made any easier through this formula, as computing 
the functions is difficult. However, this illustrates the general fact that one does not always need to 
compute the generating function for lattice point counts in every chamber to compute Hg{z)- Specifically, 
when deriving Proposition 6.1, we used the fact that the automorphism group of Ki^^-i is &k-i and the 
automorphism action of the group permutes the regions of by linear transformations. Thus, we were 
able to compute the open Ehrhart series for by calculating the Ehrhart series for one region for each 
orbit of the group action and sum these series with weights according to the size of each orbit to form 
the open Ehrhart series for ^jt- When we consider the inside-out polytope := ([0, \]'^,J^^{L)) for the 
Laplacian L of an arbittary graph G with automorphism group Aut(G), the automorphism group will act 
on the regions and the open Ehrhart polynomial for will have a decomposition as a weighted sum of 
Ehrhart polynomials for regions representing orbits of the actions. 

Our final goal in this section is to study the vertices of some of the ^/'s, thus yielding some insight into 
the denominators of Sk{z). First, it is a straightforward exercise to show that £Pl is a pyramid over [0, 1]"^^ 
and hence all the vertices of are integral. Not surprisingly, other regions of contain rational vertices. 
As we are about to see, for the case of ^| the presence of vertex denominators from the set { 1 , 2, . . . , A: — 1 } 
provides some explanation for the denominators occuring in the rational forms for Sk{z), k = 3,4,5,6. 
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Proposition 6.2. For any k>3, the following are contained in the vertex set of 



k-\ 
1 














V / 
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V J 
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V J 
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Proof. It is straightforward to check that these vectors are contained in 
are actually vertices. The constraint matrices for are of the form 



1 

V 1 J 

We need only check that these 
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: 0,. — 2, let Vi be the vector 
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^ferT3i' ^' ]t=T=i' • ■ • ' fc3T3i'0> • • • ,0)^- For i > 0' the tight inequalities 
for V, are given by the first row of the above system, the third through (3 + / — 1 rows of the system, and the 
rows of the identity matrices where equality holds, i.e., 1 or entries in v, . If i = 0, the tight inequalities for v,- 
are given by the first row of the above system and the rows of the identity matrices where equality holds. The 
matrix of coefficients for the tight inequalities is thus a square kxk matrix whose determinant is equal, up to 
sign, to the / x / minor L^t,, of with rows and columns indexed by the vertices Z,- := {1,3,4, . . . ,3 + /— 1} 
when i > and {1 } when i = 0. To show that v,- is a vertex, we need only show that det(L^,;) ^ 0. 

The non-vanishing of det(L;t,() is trivial when i = 0. For the other cases, this follows from standard 
results in algebraic graph theory, specifically Lemmas 7.4 and 7.5 of [4]. To summarize the argument, the 
Cauchy-Binet theorem implies that 



(2) 



det(L,,,-)=I^[det(ai(X,-,7))]2, 



where the sum is over all subsets Y of the edges of ^i,fe-i such that | = |7| and di {Xi,Y) is the submatrix 
of the matrix corresponding to the boundary map for Ki k_i indexed by the vertices X, and edges Y. That 
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this sum is non-zero is easily checked by applying [4, Lemma 7.4]; in the notation of the lemma, to choose 
a Y that satisfies the conditions of the lemma, we have to include in Y all the edges induced by Xi and we are 
then allowed to choose one more edge to make = |y |. As we begin with a tree, this choice of Y satisfies 
the conditions of the lemma. Hence det(5i {Xi,Y)) = 1 for some Y. As (2) is a sum of squares and one term 
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